Continuous magnetic and structural phase transitions in Fei +y Te. 
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We report a sequence of continuous phase transformations in iron telluride, Fei+^Te (y ~ 0.1), 
which is observed by combining neutron diffraction, magnetic susceptibility, and specific heat mea- 
surements on single crystal samples. While a gradual increase of magnetic scattering near the 
wave vector (0.5, 0, 0.5) is seen below T ~ 70 K, a temperature where the discontinuous first order 
magneto-structural phase transition is found in systems with small y (< 0.06), the reduction of the 
lattice symmetry in Fei.iTe only occurs at T s « 63 K. Below Tjv ~ 57.5 K the long-range mag- 
netic order develops, whose incommensurate wave vector Q m varies with temperature. Finally, at 
T m < 45 K the system enters the low-T phase, where Q m is locked at ~ (0.48,0,0.5). We conclude 
that these instabilities are weak compared to the strength of the underlying interactions, and we 
suggest that the impact of the Fe interstitials on the transitions can be treated with random-field 
models. 

PACS numbers: 71.27. +a 74.20.Mn 74.70.Xa 75.25.-j 78.70.Nx 



I. INTRODUCTION 

Iron telluride is an end member of the simplest chalco- 
genide family of iron-based high-temperature supercon- 
ductors (HTSC), Fei+yTei-zSea;. It becomes supercon- 
ducting upon partial (or full) isoelectronic substitution 
of Te by Se»i~— Although the highest critical tempera- 
ture for FeTei_ x Se x is only T c « 14.5 K, it increases 
to above 30 K in K K Fe2Se2, or under pressured The 
FeTe crystal structure consists of a continuous stacking 
of square-lattice layers of iron atoms, each sandwiched 
between the two half-density layers of bonding chalcogen 
atoms, which is the basic structural motif for all iron- 
based superconductors. The Te atoms, which tetrahc- 
drally coordinate the Fe sites, occupy alternate checker- 
board positions above and below the Fe layer, so that 
the resulting unit cell contains two formula units. In this 
quasi-two-dimensional structure, FeTe layers are held to- 
gether only by weak Van der Waals forces. Crystallo- 
graphic stability is improved if some amount of extra 
Fe atoms is incorporated between the layers, which has 
important consequences for the low-temperature phases 
observed in the Fei +y Te series, 0.02 < y < 0.18£-ii 

Similarly to the cuprate and the ferro-pnictide HTSC 
families, the end member Fei+ y Te has a magnetically 
ordered ground state and undergoes a structural dis- 
tortion which lowers the high-temperature tetragonal 
(HTT) lattice symmetry^ - — The physics behind these 
low-temperature phases and their relation to the super- 
conductivity are of great interest and have been the sub- 
jects of intense studyJ^— Two general trends of the 
phase diagrams were established: (i) unless there is a first 
order magneto-structural transition, the lattice distor- 
tion usually occurs at a higher temperature (T s ) than the 
magnetic ordering (Tjy), T s > Tjv, and (ii) both T s and 



Tjv are reduced upon chemical substitution, so that both 
orders tend to disappear as the superconducting state de- 
velops. While this observation suggests a close connec- 
tion between the magnetic ordering, the lattice distortion 
(LD), and the superconductivity, other studies indicate 
that there might be no causal relationship between these 
phenomena.— ~— In that case, they are simply different 
manifestations of the same complex physics underlying 
the electronic behavior of HTSC materials. The mecha- 
nisms by which composition affects the nature of these 
magnetic and structural transitions (first or second order, 
which occurs first), as well as their possible connections 
to superconductivity, are still not well understood. 

Fei +y Te provides an opportunity for investigating this 
issue. It is non-superconducting, but recent neutron 
studies have discovered that both magnetism and the 
low-T crystal structure are extremely sensitive to the Fe 
stoichiometryi&i We must note that the crucial issue of 
controlling the Fe stoichiometry, y, is also a tedious one. 
In particular, we measured the iron content, y, in sev- 
eral representative small crystals using the inductively- 
coupled plasma (ICP) method, and then performed sup- 
plementary neutron powder diffraction (NPD) measure- 
ments on specimens obtained by grinding single-crystal 
pieces^ We found that y obtained by refining the occu- 
pancy of the interstitial site in NPD is typically about 3% 
higher than the chemical Fe content obtained with ICP, 
y = jjicp ~ Vnpd — 0.03, for y > 0.04. Similar discrep- 
ancy was also observed by other groups^ It might imply 
that a certain amount of Fe vacancies exist, along with 
the interstitials. Here, we use the chemical Fe content 
as measured by the ICP as the appropriate notation for 
y in the Fei +2/ Te formula. This has to be kept in mind 
when comparing our results with other studies, such as in 
Refs 0-0, where different determinations of y have been 
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used. We keep their original notations when discussing 
these results. Finally, we have also found that different 
crystals from the same growth can have different y, thus 
calling for extreme caution when preparing powder spec- 
imens, as well as for the need of properly characterizing 
y for each sample. 

At low y (< 0.06), Fei +y Te undergoes a first-order 
magneto-structural phase transition, where the HTT 
PA/nmm lattice symmetry is reduced to the mono- 
clinic P21/m, and a "bi-collinear" magnetic order ap- 
pears with propagation vector (1/2,0, 1/2) in PA/nmm 
reciprocal lattice units'— The main peculiarity of the 
"bi-collinear" magnetic ordering in Fei+yTe materials 
compared to the simple bipartite antiferromagnetism of 
other HTSC families is that it does not agree with the 
Fermi surface nesting of itinerant electrons, which cor- 
responds to the (1/2, 1/2) position in the a6-plane^r— 
Hence, band structure calculations, which tend to ac- 
count for the Fermi surface^ - — encounter difficulties 
when confronted with the broad range of experimental 
observations ^ 10 ' 32 ! 37 " — . On the other hand, the ground- 
state ordered moment, (fj,) < 2/iBr 1 ^^- although larger 
than in parent ferropnictides, is significantly less than 
4/i_B (p-B = Bohr magneton) expected in the ionic local- 
spin picture for Fe 2+ (S = 2) in the Hund's state,- 4 While 
some theories postulate that Fe could be in the non- 
Hund, intermediate-spin S=l stated such an assump- 
tion falls way short of accounting for the paramagnetic 
moment fJ> e ff ~ 4/xb implicated in the Curie- Weiss be- 
havior above 100 K 3 - 7 - 3 - 8 . Thus, in Fei +a Te we can explore 
the interaction of weak magnetic order with lattice and 
orbital degrees of freedom in a representative structure of 
Fe-based superconductors, disentangled from the effects 
of Fermi surface nesting and superconductivity. 

The phase diagram for y > 0.06 is still controver- 
sial. Bao et aZ.— who first discovered the orthorhombic 
Pnmm phase with an incommensurate magnetic order 
for y w 0.14, initially suggested that the incommemn- 
surability varies linearly with y, and the transition is 
first order. A more recent NPD study 7 - suggested that 
the first-order transition to a P21/m, "bi-collinear" com- 
mensurate phase survives until y « 0.12, at which point a 
mixed phase is observed. At higher y, the low-T phase is 
orthorhombic, with coexisting long-range helimagnetism 
and short-range spin-density-wave (SDW) order. 

Here, we report studies of well-characterized single 
crystals of Fei+ y Te, y = 0.10(1), with a variety of tech- 
niques, which establish a sequence of continuous phase 
transitions. It starts with a structural distortion at 
T s = 63(1)K, which is followed by slightly incommen- 
surate magnetic order at Tjv = 57.5(5) K. This im- 
plies a multicritical point, y c on the (y, T) phase dia- 
gram of Fei +J ,Te, with 0.06 < y c < 0.1, where the first- 
order magneto-structural phase transition turns into a 
sequence of continuous ones. While according to Ref. 7 
our samples belong to the grey "mixed" phase, we find 
no evidence for the mixed character. Instead, we find a 
well-defined sequence of phases as a function of temper- 



ature. This leads us to favor the idea that samples stud- 
ied in Ref. [7J were mixtures of different stoichiometrics, 
y, among other possibilities suggested by the authors of 
that study. This conclusion is further supported by our 
more recent additional NPD and bulk measurements of 
the (y, T) phase diagram of Fei +y Te..— They suggest that 
multicritical point at y c is indeed an intrinsic property of 
this phase diagram. Hence, our findings provide a di- 
rect connection between the composition phase diagram 
of Fei+ y Te and that of BaFe2As2 (122) derived materi- 
als, where a similar multicritical point has recently been 
focus of considerable work. 18-23 



II. EXPERIMENTAL PROCEDURE 

The Fei.iTe crystals were grown by the horizontal 
Bridgman method. 3 The crystal used for the neutron- 
scattering study had a mass m = 18.45 g and a mosaic 
of 2.2° full-width at half-maximum (FWHM). 

Our neutron measurements were carried out using the 
ARCS direct-geometry, time-of-flight (TOF) spectrom- 
eter at the Spallation Neutron Source, Oak Ridge Na- 
tional Laboratory. The instrument was operated in the 
Laue-diffraction-type mode^ where a quasi-white neu- 
tron beam with a broad band of incident neutron en- 
ergies centered around Ei m 300 meV was selected by 
the pre-monochromating Tq chopper. The Fei.iTe crys- 
tal was mounted on an aluminum holder attached to the 
cold head of the closed-cycle refrigerator in the ARCS 
evacuated scattering chamber. The crystal's c-axis was 
aligned in the horizontal plane, at w 45° to the incident 
neutron beam, while the a-axis was at about 24° to the 
horizontal plane. The detector signal is dominated by 
the elastic processes (diffraction), where the scattering 
angle is determined by the incident neutron energy (or 
wavelength, Aj) and the rf-spacing of the set of crystal 
planes involved in reflection, in accordance with Bragg's 
law, Xi = 2d sin 9. Such a measurement is particularly 
well suited for studying the relative temperature evo- 
lution of structural and magnetic scattering, which are 
both present in the diffraction pattern at each T. 

Bulk magnetization, M, and heat capacity, C p , were 
measured on several single crystals with masses from 
7.5 mg to 26 mg, cleaved from the same growth boule 
as the large crystal used for neutron studies. C p was 
measured using the relaxation method implemented in 
the Physical Properties Measurement System (PPMS) 
by Quantum Design (QD). Crystals were attached to a 
silver sample holder using apiezon grease, whose contri- 
bution to C p was measured and subtracted. Its uncer- 
tainty was the main source of the dominant systematic 
error. The non-magnetic lattice phonon specific heat was 
estimated using two different algorithms. First, we used 
the C P (T) of non-magnetic ZnTe,— with temperature re- 
scaled by the ratio of the effective Debye temperatures, 
6d, of Fei.iTe and ZnTe, obtained by fitting their C P (T) 
in the T > 8d range to a single Debye function (short- 
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Figure 1: (a) Temperature dependence of the ZFC static mag- 
netic susceptibility of the Fei.ogTe single crystal measured in 
DC magnetic field B — 0.1 T applied in the ab— plane, Xab[T), 
and along the c— axis, Xc(T). (b) inverse susceptibility, illus- 
trating the Curie- Weiss linear asymptotics at high temper- 
atures. Dashed lines are Curie- Weiss fits, (c), (d) The ZFC 
and the field-cooled (FC) susceptibilities and their derivatives. 
The shaded region shows phase with temperature-dependent 
magnetic incommensurability Q m , the hysteretic range of the 
lock- in transition is light-shaded. 

dashed lines in Fig. [2]). Secondly, we used an equal- 
weight sum of two Debye and one Einstein functions, as 
described in Ref. H|| for the FeSeo.sTeo.s case, re-scaled 
in a similar way (solid lines in Fig. [2]). 

The static magnetic susceptibility, \ — M/H, was ob- 
tained from M measured using a QD Magnetic Properties 
Measurement System, in a dc magnetic field (j,qH = 0.1 T 
applied in the afr-plane, Xabi or along the c-axis, \c- 
The iron content, y — 0.10(1), was measured in these 
and several other representative small crystals, using 
the inductively-coupled plasma method. The 10% error 
bar on y results from the scatter of different ICP mea- 
surements. We performed supplementary neutron pow- 
der diffraction measurements on specimens obtained by 
grinding single-crystal pieces 2 ! and established that y ob- 
tained by refining the occupancy of the interstitial site is 
about 3% higher than the chemical Fe content obtained 
with ICP, y ee yi C p « y NPD - 0.03. 



III. RESULTS AND DISCUSSION 

Figure [Ija) shows the temperature dependence of the 
magnetic susceptibility measured upon heating the zero- 
field-cooled (ZFC) sample. For T > 100 K, it obeys the 
Curie- Weiss (CW) law, 



Xa (T) = N A 



3k B (T- 9 



CW,a) 



(1) 



where Na is Avogadro's number, kg is Boltzmann's con- 
stant, a = ab, c indexes the field direction, and /z® 
and 6cw, a are the effective paramagnetic moment and 
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Figure 2: (a) ZFC (light-filled) and FC (dark-filled symbols) 
specific heat capacity of Fei.ogTe. The solid curve is the esti- 
mated lattice vibrational contribution, the dashed curve is the 
Debye fit in the T > 130 K range, the inset shows the corre- 
sponding net magnetic specific heat, (b) The total (light-filled 
symbols) and the net magnetic (connected points with error 
bars) entropy. Curves show the estimated lattice contribu- 
tion. Inset expands the region near Tn = 57.5(5) K, which 
is marked by the arrows, (c) and (d) show the net magnetic 
heat capacity and the magnetic entropy, respectively. Sym- 
bols and error bars connected by the dashed lines result from 



the Curie- Weiss temperature. This behavior is consis- 
tent with previous studies ^ 37 ' 38 and is best revealed by 
plotting 1/Xab,c as m Fig. 1(b). CW fits to our data 
in the range @cw T < 300 K yield large effective 
magnetic moments, /i® ff = g a ^B\J ' S(S + 1), which are 
consistent with Fe atoms having local spins S — 3/2 with 
slightly anisotropic Lande g— factors, g a b ~ 2.15 and g c ss 
2.02. The CW temperatures are negative, correspond- 
ing to a dominant antiferromagnetic interaction, and are 
also very slightly anisotropic, with 0cw, ab — —165 K 
and Ocw,c = —157 K. Remarkably, the magnitude of 
&B@cw,a is much smaller than the bandwidth of mag- 
netic excitation o 42 ' 43 , which indicates strong competition 
between ferro- and antiferromagetism. The low-T suscep- 
tibility and its derivative in Figure [He), (d) clearly reveal 
a continuous magnetic phase transition at Tn — 57.5(5) 
K. x(T) has a cusp at Tn, rather than a first-order- 
like discontinuity, which is further corroborated by the 
A-like (step) singularity in d\jdT. It is followed by an- 
other magnetic transition, at a lower temperature, show- 
ing ZFC-FC hysteresis in the 35-45 K range. 

The continuous, second-order nature of magnetic or- 
dering at Tn is further corroborated by the heat capacity 
in Fig. [21 which shows a A— type singularity. The change 
in magnetic entropy associated with the long-range mag- 
netic ordering at Tn is very small, Fig. (2]^b),(d). It does 
not even reach kg In 2 per Fe, which would correspond to 
freezing of a single Ising degree of freedom per Fe. This, 
together with a rather small ordered magnetic moment 
observed by neutron diffraction, (/i) < 2[iBr^ , ]10]43 
shows that long-range order (LRO) is weak. The most 
likely reason for the weakness of magnetic order is frus- 
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tration arising from the competition of ferro- and antifer- 
romagnetic, nearest- and further-neighbor interactions, 
which is also indicated by the small 0cw ■ The difference 
between the ZFC and the FC C P (T), as well as that mea- 
sured upon heating the field-cooled sample in a magnetic 
field of B = 7 T, is very small, suggesting that the hys- 
teretic transition observed in susceptibility at 35-45 K is 
some sort of spin realignment, involving negligible change 
in magnetic entropy. 

Figure [3J presents an overview of our neutron diffrac- 
tion data at two temperatures, T = 9K and 80 K, in the 
form of Laue patterns on the detector bank. Each de- 
tector pixel is parameterized by the unit vector n(0, </>) 
specifying its direction from the sample position. Al- 
though rather un-intuitive alignment of the sample recip- 
rocal space with respect to the incident beam direction 
(ps 24° rotation around plus ps 45° rotation around the 
vertical axis) results in an un-obvious pattern of equiv- 
alent Bragg peaks, reflections can be easily identified by 
their d— spacings. Panels (a), (b) show intensities corre- 
sponding to elastic scattering for the range of d— spacings 
from 1.75 A to 1.9 A, which is dominated by the lattice 
nuclear scattering. As magnified in the insets, a single 
Bragg spot corresponding to the (2, 0, 1) lattice reflection 
seen at 80 K is clearly split into two spots at T = 9 K. 

Intensities in the bottom panels (c), (d) of Fig. [31 cor- 
responding to 5.75 A< d < 7.75 A, are mainly magnetic 
scattering. At the base temperature of 9 K it appears 
in the form of well-defined peaks near (±1/2, 0, ±1/2) 
positions. Their significant angular size on detector re- 
sults from the combination of deteriorating experimental 
resolution at small wave vectors (large d— spacings) and, 
to some extent, from contribution of quasielastic and in- 
elastic diffuse scattering, which is also collected in this 
measurement. At T — 80 K only a small intensity modu- 
lation highlighting former peak positions remains, while 
magnetic intensity appears distributed over a ring of scat- 
tering. Such pattern is characteristic of a liquid. 

The competition between different magnetically or- 
dered states, confirming the frustrated nature of mag- 
netic interactions, is further revealed by the temperature 
dependencies of quasi-Laue neutron data in Figs. U and 
[5j They show intensities, obtained by integrating pat- 
terns, such as in Fig. [3J over the whole angular range 
of the detector array, as a function of d— spacing. Mag- 
netic peak profiles near the Q = (1/2,0,1/2) position, 
Fig. Bib) and Fig. [SJc,d), show complex structure evolv- 
ing with temperature. Broad diffuse scattering around 
d[Q ps (0.4,0,1/2)] ps 7.6 A and a narrower peak near 
d[Q w (0.51, 0, 1/2)] ps 6.4 A coexist for T > 70 K. Below 
pa 70 K they yield to a peak at d[Q ps (0.485, 0, 1/2)] ps 
6.65 A, whose intensity grows roughly linearly with the 
decreasing T. Then, below ps 57.5 K, a new dominant 
peak emerges at d[Q ps (0.46,0, 1/2)] ps 6.9 A. Its posi- 
tion changes upon cooling, and finally saturates at d[Q ps 
(0.48,0, 1/2)] ps 6.7 A below ps 45 K. At the same time, 
a smaller peak appears at d[Q ps (0.52,0, 1/2)] pa 6.3 A. 

While very little, if any of these behaviors could be 
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Figure 3: Contour maps of the scattered neutron intensity on 
the detector bank, as a function of the scattering direction, 
n, at T = 9 K (a), (c) and T = 80 K (b), (d). (n x ,n y ) are 
projections of the unit vector n pointing from sample to a de- 
tector pixel onto the plane perpendicular to the incident neu- 
tron beam direction, ki. Top panels, (a), (b), show the range 
of d— spacing dominated by structural scattering, and reveal 
the splitting of the (201) nuclear Bragg reflection on cooling. 
Bottom panels, (c), (d), show the large— d range dominated 
by magnetic scattering, which transforms from a set of well- 
defined peaks at 9 K, to a ring-like feature typical of a liquid 
at 80 K. Insets show parts of the same data (marked by ar- 
row) in angular coordinates, (6>,<jf>), where n x — sin 9 cos <j>, 
n y — sin(9sin</>. The intensities have been window- averaged 
over the window of 3 x 3 detector pixels. 



identified in the integrated magnetic intensity shown in 
Fig. [3Jc) , they are clearly observable in the temperature 
dependence of the peak maximum intensity, Fig. EJa), 
and its position and width, Fig. [4f d) , (e). In particu- 
lar, the appearance of the new dominant magnetic com- 
ponent at ps 57.5 K is most clear from an abrupt shift 
of the maximum intensity position in Fig. HJd). Since 
it appears at the temperature where magnetic order is 
observed in susceptibility and heat capacity data, we 
identify this peak as a magnetic Bragg reflection asso- 
ciated with magnetic LRO. Fitting the peak intensity for 
T > 50 K in Fig. (5^a) to an order-parameter-like de- 
pendence, I(T) ~ (T N - T) 2/5 , we obtain T N = 57.5(5) 
K and 2/3 = 1.0(1), consistent with the mean-field, lin- 
ear I(T) behavior. The average peak position and its 
effective width in Fig. H^d) and (e) are governed by the 
structured, multi-component nature of magnetic scatter- 
ing and have very large error bars at high T, where the 
net magnetic intensity is small. Note that error bars on 
the maximum intensity position, which are shown by the 
closed symbols in Fig- Hl^cL) , are much smaller, of the or- 
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Figure 4: (a) Representative scans through the lattice Bragg 
peak (2, 0, 1) and (b) the magnetic scattering near Q = 
(0.5,0,0.5). The data at different T are vertically offset for 
presentation, (c) The integral intensity of the magnetic (filled 
blue symbols) and the lattice (open black symbols) scatter- 
ing; (d) the maximum intensity position (filled blue symbols) 
and the intensity-weighted average position (open blue sym- 
bols) of magnetic peak and the intensity-weighted average of 
the lattice (2,0, 1) peak (open black symbols); (e) the width 
(intensity-weighted mean-square deviation of the measured 
points from the peak center) of magnetic (filled blue symbols) 
and the lattice (open black symbols) scattering. 
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Figure 5: (a) The (maximum) peak intensity of the lattice 
(open) and the magnetic (filled symbols) scattering, (b) The 
integral intensity (filled symbols) and the width (intensity- 
weighted mean-square deviation of the measured points from 
the peak center, open symbols) of the (2, 0, 1) lattice Bragg 
peak. Solid lines are fits to ~ \T s ,n — T| 2/3 dependence. Dot- 
ted vertical lines show T s and Tjv. (c) and (d) show rep- 
resentative magnetic scattering data at several temperatures 
on linear and logarithmic scale, respectively. Dashed vertical 
lines show d— spacing for Q = (0.5, 0, 0.5). 

der of the symbol size. 

Since fitting the overlapping peaks is sensitive to fit- 
ting ranges and constraints, peak parameters in Figs. 0] 



and [3] were evaluated directly from the measured inten- 
sities. The peak integral intensity, position, and width, 
were obtained by numerical integration of the measured 
intensity, the intensity-weighted position, and the mean- 
square deviation, respectively, upon subtracting the lin- 
ear background interpolated between the edges of data 
ranges shown in Fig. 0Ja) and (b). Parameters of the 
(2, 0, 1) structural Bragg reflection, whose splitting at low 
temperature reveals the lattice distortion that reduces 
the symmetry from HTT, Fig. HJa), were obtained in the 
same way. The position of magnetic and lattice peaks 
in Fig. [3](d) are given relative to the nominal positions, 
d[Q = (1/2,0, 1/2)] and d[Q = (2,0, 1)], respectively, in 
the HTT lattice with a = 3.813 A, c = 6.24 A. 

While the evolution of the (2,0, 1) reflection between 
a single peak at 100 K and a two-peak structure at « 9 
K seems gradual (mainly due to the experimental reso- 
lution), both the integral and maximum peak intensities 
in Fig. |4jc) and Fig. [5ja), as well as the peak width in 
Fig. [5jb) , immediately reveal the structural phase tran- 
sition at w 63 K. The width directly probes the order 
parameter - the splitting of the (2,0,1) peak. Fitting 
it to an order-parameter-like dependence in T > 48 K 
range, we obtain T s — 63(1) K and 2/3 = 1.0(1), again 
consistent with the mean-field behavior. The intensity 
is an indirect probe, and its temperature dependence is 
governed by the combination of splitting and extinction, 
which is why the two dependencies in Fig. [SJb) differ. 
The lattice Bragg intensity measured on a large single 
crystal is reduced as a result of neutron beam extinction 
within the crystal. It is therefore very sensitive to T s , 
where the crystal's mosaic structure changes due to the 
appearance of domains associated with the lowering of 
the HTT symmetry ? ?. While our present data do not 
allow us to distinguish between the monoclinic P2\/m 
and the orthorhombic Pnmm structures^ our supple- 
mentary powder diffraction measurements on a sample 
with similar y indicate the P21/m structure. 



IV. SUMMARY AND CONCLUSIONS 

To summarize, our analysis of neutron Bragg scat- 
tering establishes the structural phase transition, low- 
ering the HTT lattice symmetry at T s — 63(1) K, as 
the first instability which occurs in Fei.iTe upon cooling. 
This transition is continuous and involves only very small 
structural changes: the (2,0, 1) peak splitting, Sd/do, at 
9 K is only ~ 0.5 %. Moreover, this transition does 
not show up neither in heat capacity, Fig. [21 nor in bulk 
magnetic susceptibility, Fig. [1] indicating negligible en- 
tropy change and nearly complete decoupling between 
this structural change and the magnetic order. Close ex- 
amination of the magnetic scattering in Figs. [5] and [S] 
further corroborates this observation: its complex tem- 
perature evolution shows no visible anomaly at T s . 

Thus, the lattice distortion, which occurs in Fei.iTe at 
T s = 63(1) K, is not induced by the long-range antifer- 



6 



romagnetic ordering. Neither does it immediately lead 
to magnetic LRO, which only follows at a « 10% lower 
temperature, Tjv = 57.5(5) K. Magnetic order is weak. 
It accounts for freezing of only < 25% of (fcsln4/Fe) 
paramagnetic entropy of S = 3/2 spins, implicated in 
the CW behavior of magnetic susceptibility at T > 100 
K. While the present quasi-Laue data is not suitable 
for the absolute normalization and refinement of the 
LRO moment involved in Bragg scattering, complemen- 
tary monochromatic beam measurements 4 - indicate an 
ordered moment ~ 1.4/ig in our Fei.iTe crystal. This 
agrees with rather small values of the ordered magnetic 
moment, (/x) < 2/ig, observed by neutron diffraction in 
Fei +! ,Te material o 6 i 7 i 9 i 10 . 

Magnetic LRO in our sample is not simple "bicollinear" 
type, but is an incommensurate structure, whose pe- 
riod varies with temperature. Upon cooling, it un- 
dergoes what looks like a lock-in transition with sig- 
nificant FC-ZFC hysteresis in magnetic susceptibility, 
as its propagation vector saturates at a low-T value 
Q rs (0.48,0,1/2). In real space, such structure can 
be visualized by introducing a certain amount of "tri- 
collinear" defects, whose density decreases with temper- 
ature. These defects can be viewed as randomly in- 
serted lines of corner-sharing square plaquettes with fer- 
romagnetically co-aligned spins. It was recently found 
that such plaquettes govern low-energy spin dynamics in 
Fei.iTe 43 . Such proliferat ion of defects highlights mag- 
netic frustration, where a number of different magnetic 
ground states, including the plaquette state, have nearly 
equal energy, 49 This situation is further demonstrated by 
the temperature evolution of the magnetic neutron in- 
tensity, where, upon cooling, one observes competition 
of a number of states with different propagation vec- 
tors. We note that incommensurate magnetic orderings 
were also observed in a number of Se-doped compositions, 
Fe 1+!/ Tei_ a Se x ,29-ii 

Finally, our observations suggest that the weak lattice 
distortion and antiferromagnetic LRO are only symp- 
tomatic of the much stronger interactions driving the 
low-temperature physics in iron telluride. This leads us 
to question the applicability of traditional theoretical ap- 
proaches, which treat electronic and magneto-structural 
properties in terms of expansions around the low-T or- 
dered states. The low-temperature physics in Fei + j,Te, 
and, perhaps, in ferropnictides, too, is likely governed 
by rather high-energy degrees of freedom, such as the 
temperature-dependent orbital hybridization and inter- 
action of local spins with itinerant electrons , i 50 Upon 
being properly integrated, these interactions should ren- 
der an effective low-energy model governing the magneto- 
structural transition. Is there such a model, which could 
explain the extreme sensitivity to small amount of dop- 
ing and the multi-critical phase diagrams common to 
Fei +l) Te and 122 pnictides? 

The effect of the small Fe off-stoichiometry, y, on the 
electronic structure of Fei+ y Te being negligible^, Fe in- 
terstitials/vacancies can simply be viewed as a source 



of random magnetic and electric fields, which frustrate 
magnetic correlations, locally impact lattice distortion, 
and change the inter-layer structural and magnetic cou- 
plings. An effective description of the spontaneous HTT 
symmetry breaking in layered systems under these con- 
ditions is provided by the quasi two-dimensional (2D) 
anisotropic random field Ising model (ARFIM) consid- 
ered by Zachar and Zaliznyak^i While this model was 
originally proposed for the superlattice formation associ- 
ated with charge ordering in layered perovskites at half- 
doping, it is straightforwardly transplanted to the case 
of an orthorhombic distortion of the HTT lattice. The 
two possible choices of distortion at a given lattice site, 
which are related by the 90° rotation, now play the role 
of an Ising variable. The effective ferromagnetic interac- 
tion accounts for the strain energy, arising where these 
two different states are adjacent to each other, and is 
strongly anisotropic (quasi-2D) in layered systems. The 
most important property of the ARFIM model is in its 
extreme sensitivity to very small amounts of disorder in 
the quasi-two-dimensional limit, which is inherited from 
the disordered nature of RFIM in 2D^i In fact, recent 
studies have found that ARFIM can account for signa- 
tures of nematicity observed in cuprates^ and thus is 
probably also a good starting point for understanding 
similar findings in ferropnictides 15 . 

A natural generalization of the ARFIM 51 to cases 
where the monoclinic lattice distortion and/or strong 
magneto-elastic coupling result in four choices of the 
ground state, such as for the magneto-structural tran- 
sitions in 122 ferropnictides and 11 chalcogenides, is pro- 
vided by the four-state (q = 4) anisotropic random field 
Potts model (ARFPM). Similarly to the ARFIM, the or- 
dering in the ARFPM is governed by the 3D fixed point 
of the random field Potts model (RFPM). In RFPM, 
however, it is of the first order for small disorder^ A 
multicritical point, where the transition changes from 
first-order into a continuous one with the increasing ran- 
dom field strength, is expected in this case^ Hence, such 
a model could provide natural explanation for our cur- 
rent results, suggesting that such a multicritical point 
exists in Fei +y Te, as well as for the recent findings in 
122 ferropnictides ^r21 
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